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We report combined soft and hard x-ray scattering studies of the electronic and lattice modulations 
associated with stripe order in Lai.875Bao.i25Cu04 and Lai.4sNdo.4Sro.i2Cu04. We find that the 
amplitude of both the electronic modulation of the hole density and the strain modulation of the 
lattice is significantly larger in Lai.87sBao.i25Cu04 than in Lai.4gNdo.4Sro.i2Cu04 and is also better 
correlated. The in-plane correlation lengths are isotropic in each case; for Lai.875Bao.i2sCu04, 
£hoie = 255±5 A whereas for Lai.4 8 Ndo.4Sro.i 2 Cu04, £ hole = 111±7 A. We find that the modulations 
are temperature independent in Lai.875Bao.i2sCu04 in the low temperature tetragonal phase. In 
contrast, in Lai.4gNdo.4Sro.i2Cu04, the amplitude grows smoothly from zero, beginning 13 K below 
the LTT phase transition. We speculate that the reduced average tilt angle in Lai.875Bao.i25Cu04 
results in reduced charge localization and incoherent pinning, leading to the longer correlation length 
and enhanced periodic modulation amplitude. 



PACS numbers: 71.30.+h, 61.10.Nz, 71.27.+a 

One dimensional modulations of the lattice and of the 
charge and spin density appear to be both ubiquitous in 
the cuprates and intimately tied up with the physics of 
these materials. For example, superconducting phase co- 
herence is suppressed in systems in which charge and spin 
stripes are statiJ^' while dynamic stripes may promote 
pairing^!. Further, stripe order coexists with 2D super- 
conducting correlations above the bulk superconducting 
phase^l and indeed it has been argued that broken rota- 
tional symmetry may be a defining feature of the pseudo- 
gap phase - so-called "nematic" electronic order, a view- 
point for which there is some evidence^ . Understanding 
these modulations is therefore a prerequisite to gaining a 
detailed knowledge of the physics of the cuprates. 

In their full incarnation, these density waves - "stripes" 
- consist of three concomitant modulations; of the spin 
density, of the charge density and of th e latt ice it- 
self. In a few cases, namely La2- :r Ba a ;CuO i^ 7 -^ ^ and 
(La,R) 2 - 2; Sr x Cu04, where R=Nd 13 -, or ErP^ , these 
modulations are stabilized by the symmetry of the low 
temperature crystal structure to form static structures. 
It is here, where the respective order parameters are well 
developed, that one can best characterize them. The 
basics of the static structures are known: they are com- 
prised of hole-rich anti-phase antiferromagnetic domain 
walls running along the Cu-0 bond direction. In real 
space the spin density then has twice the wavelength of 
the charge and lattice modulations, and in a scattering 
experiment, these give rise to satellites at (0.5 ± 5,0.5,0) 
and (±25,0,0.5), and symmetry related positions, respec- 
tively. However, basic questions such as the relative am- 
plitudes in the different compounds and the coupling be- 



tween these modulations remain unanswered. 

Here we report a combined soft and hard x-ray scat- 
tering study of the low temperature stripe order in 
Lai.875Bao.i25Cu04 (LBCO) and Lai.4sNdo.4Sro.i2Cu04 
(LNSCO). In particular, soft x-ray scattering at the pre- 
edge of the O K-edge (529.3 eV) probes the spatial vari- 
ation of the oxygen 2p partial density of states at the 
Fermi leve l, which is related to the mobile charge car- 
rier density^- * 16 * 17 !, while hard x-ray (11.3 keV) scatter- 
ing probes the associated lattice distortion. We find that 
the stripe correlation lengths are isotropic in the Cu-0 
plane in both systems, but LBCO is a factor of two better 
correlated. Further, by careful normalization to struc- 
tural Bragg peaks we are able to quantify the relative 
size of the electronic modulation and the lattice distor- 
tion in these two cases. We find that the amplitude of 
the electronic modulation is ~ 10 times larger in LBCO 
than in LNSCO, while the lattice distortion is only ~ 4 
times larger. We further find that the modulations are 
temperature independent in LBCO once the LTT phase 
is entered, suggesting that perhaps the intrinsic order- 
ing temperature is significantly higher. In contrast, in 
LNSCO, the amplitude grows smoothly from zero, be- 
ginning 13 K below the LTT phase transition. 

Samples of both LBCO and LNSCO were grown using 
the floating zone method, at Brookhaven National Labo- 
ratory, USA and the University of Tokyo, Japan respec- 
tively. Both samples proved to be of high crystal quality, 
with mosaic spreads of < 0.02° measured at 11.3 keV. 
Throughout this paper both systems are indexed with 
the high temperature tetragonal (HTT) (IJ^/mmm) unit 
cell. 
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Soft x-ray diffraction experiments were carried out on 
the XI A2 beamline at the NSLS, Brookhaven National 
Laboratory, USA using a six-circle in-vacuum diffrac- 
tometer. The samples were cleaved ex-situ to reveal sur- 
faces with a [001] htt surface normal and mounted such 
that the [001]htt and [100]htt directions lay in the scat- 
tering plane. Experiments were performed in a vertical 
scattering geometry with er-incident x-rays, i.e. the elec- 
tric field of the incident x-ray photons was always along 
the HTT 6-axis, within the CuC>2 planes. The samples 
were cooled in a He flow cryostat to a base temperature of 
12 K. Scattered x-rays were detected with an in-vacuum 
CCD. During the measurements the beam footprint was 
smaller than the sample. Data were collected by per- 
forming a — 28 scan, collecting a single CCD image at 
each datum. Each CCD image was then converted to 
reciprocal space coordinates and binned onto a regular 
grid^. Part of the soft x-ray experiments were performed 
at BESSY with a diffractometer described in Ref. UM 

Hard x-ray diffraction experiments on the same sam- 
ples were conducted on beamline X22C at the NSLS. The 
incident photon energy was 11.3 keV and scattered x- 
rays were detected by means of a point-detector, with the 
resolution defined by slits. Experiments were conducted 
in a vertical scattering geometry, with the [001] htt and 
[100] htt directions lying in the scattering plane. For 
both LBCO and LNSCO the instrumental resolution was 
identical. 

We first discuss the characterization of the electronic 
modulation, beginning with LBCO. In Fig. [I] we show 
cuts through reciprocal space in the HK- (bottom) and 
HL- planes (top) for LBCO. These data were taken with 
Ei=529.3 eV, an energy corresponding to the Mobile 
Carrier Peak (MCP) in the oxygen K-edge absorption^. 
At this energy the scattering from the charge stripes is 
maximized™ We extract a correlation length from these 
data by taking line cuts through the three dimensional 
grid and fitting to Lorentzian-squared lineshapes. We 
find that the stripe correlations within the a-b plane are 
isotropic with = 255 ± 5 A, where £ = HV ^ HM ■ Mea- 
surements at the Cu L3 edge revealed a similar correla- 
tion length. We note that because of the superior resolu- 
tion of the present set up, this represents the highest pre- 
cision measurement to date of this correlation length. It 
is consistent with previous work which found £ ~ 200 A 
at this energy^. As has been previously observed, the 
stripes are weakly correlated in the c-axis direction, giv- 
ing rise to the uniform streak of intensity along L in 
Fig. b 

Figure [2] shows a comparison in the HK- plane, along 
with line cuts along the H direction, between LBCO and 
LNSCO. Immediately apparent is the far poorer correla- 
tion length in LNSCO, where we find that the correlation 
length is again isotropic but in this case £# = 111±5A. 

In both cases, the wavevector in the H direction, is 
0.24 r.l.u. and is therefore incommensurate. These mea- 
surements were performed using an orientation matrix 
defined in-situ by the (002)htt and (IOI)htt Bragg re- 
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FIG. 1. Soft x-ray resonant scattering intensity in the vicinity 
of the (0.24, 0, L)htt stripe order peak for LBCO in the a) 
(H, 0,L)htt plane and b) (H, K, 0.624)htt plane in recipro- 
cal space. Data were taken at T=15 K, and Ei=529.3 eV. 
At this energy, the scattered intensity reflects the modulation 
of the oxygen density of states in real space. The intensity 
map reveals that the stripe correlations are isotropic in the 
a-b plane and are two dimensional. 



flections, and as such their values are reported with high 
confidence. Previous studies have found similar incom- 
mensurabilities at these dopings (see e.g. Refs. QUI H21 
and lT9")) . We return to this point later. 

The result that the charge stripes display isotropic 
correlation lengths is perhaps surprising given that the 
stripes in a single Cu02 layer are highly anisotropic. 
While four-fold rotational symmetry is recovered when 
averaging over the whole crystal, since successive sheets 
along the c-axis are rotated by 90° with respect to each 
other, these present experiments only observe charge 
stripes with a propagation wavevector along one direc- 
tion in reciprocal space - by virtue of the diffraction con- 
dition. There is therefore no a priori reason to expect 
these correlations to be isotropic. One possible expla- 
nation is that nearest-neighbor c-axis correlations couple 
the charge stripe correlations on alternating layers. It 
has been suggested, for example, that such inter-layer 
coupling could occur between the ID metallic behavior 
of a stripe and the q-vector in the neighboring stripe^. 
Unfortunately, from our diffraction measurements it is 
not possible to say anything about nearest-neighbor c- 
axis coupling due to the 90° rotation of the stripes. A 
second possibility is that the isotropy is driven by poorly 
screened (long range) Coulomb interactions, as has been 
suggested in the context of a ladder compound^. Finally, 
and most prosaically, it is possible that the charge stripe 
correlations are limited by intrinsic crystallographic dis- 
order, for example due to defects arising from randomly 
distributed dopant ions^. As such order is random, one 
would expect it to be isotropic, and that it would there- 
fore cause an isotropic correlation length. We will return 
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to the role of crystallographic order again shortly in the 
discussion of modulation amplitudes. 

We are also able to compare the peak intensities in 
the two systems. To do so, the intensities shown in 
Fig. [2] have been normalized to the integrated intensity 
of the respective (002)htt Bragg reflections (measured 
at 1 keV) and corrected for their different structure fac- 
tors. The integrated intensity, V, of a superlattice re- 
flection is proportional to the product of the peak inten- 
sity and the widths in three orthogonal directions, V oc 
^Pcak x T# x Fk xTl- As discussed, Tk = r#, and we 
make the assumption that the correlation lengths along 
L are equal in the two materials. The ratio of integrated 
intensities is then Vjg le = (Ip ak A x r^)/(J Pea k b x T 2 b ). 
Calculating from the data presented in Table [I] we find 
that the ratio of integrated intensities, corrected for all 
experimental factors is V™j = 97±10, i.e. the electronic 
charge stripe modulation in LBCO is ~ 9.8 ± 0.5 times 
larger in amplitude than in LNSCO. 

This factor of almost fO is a surprisingly large num- 
ber. Taken at face value it implies a very large dif- 
ference in the value of the charge order parameter in 
these two systems and appears to contradict expectations 
based on measurements of the ordered magnetic moments 
in the two systems. /iSR measurements, for example, 
show similar ordered moments of 0.3 /xg with similarly 
large magnetically-ordered volume fractions in LBCO 
and LNSCCW Since the magnetic order is constrained 
by the charge order, it would seem counterintuitive for 
the charge order parameters to be so different when the 
magnetic order parameters are so similar. Further, a re- 
cent estimate of the hole concentration per Cu site in the 
charge carrier stripes in Lai.675Euo.2Sro.i2sCu04 (LE- 
SCO) resulted in a value considerably larger than 0.2 
holesP. Based on the similarity of LESCO and LNSCO, 
in terms of correlation length and temperature depen- 
dence, once might expect LESCO and LNSCO to have 
similar hole modulations. Such an expectation is then, at 
first glance, inconsistent with the present measurement, 
which suggests LBCO is 10 times larger, which would 
be an unphysical result. However, the LESCO estimate 
came from an analysis of the resonant lineshape, which 
reflects the amplitude of the modulation in the modu- 
lated regions and there are a number of reasons why a 
measurement based on the diffraction signal could lead to 
a different result. Firstly, the ordered volume of sample 
could be different for LBCO and LNSCO, with LNSCO- 
having an ordered volume up to 100 times smaller than 
that of LBCO. This seems extremely unlikely. A sec- 
ond possibility is that these diffraction measurements do 
not fully capture all the intensity of the stripe correla- 
tions and that in LNSCO there is significant intensity 
outside the detected region in reciprocal space - i.e. in 
diffuse tails. Missing this intensity would result in a much 
reduced measured hole amplitude in LNSCO compared 
with LBCO. This effect would not be reflected in reso- 
nant lineshape analysis. Such an explanation is perhaps 
consistent with the fact that the charge stripes are more 



TABLE I. Comparison of hole and lattice modulation 
wavevectors, correlation lengths and peak intensities for 
LBCO and LNSCO as measured with soft and hard x-rays 
respectively. 



Center HWHM 
[r.Lu.] [r.l.u.] xlO 
LBCO 529.3 eV 0.241 2.36 ±0.05 
LNSCO 529.3 eV 0.238 5.42 ± 0.25 
LBCO 11.3 keV- 3.37 ±0.9 

LNSCO 11.3 keV - 5.59 ± 0.9 



Peak 
" 3 Intensity 
7.82 ±0.01 
(1.54 ±0.005) x 10" 3 
0.225 ±0.002 
(4.17 ±0.18) x 10 -3 



La i 87» Ba o i25 Cu0 4 La 148 Nd 4 Sr 012 CuO 4 




FIG. 2. Comparison of scattering from stripe modula- 
tions in LBCO and LNSCO. Data were taken at T=15 K, 
and Ei=529.3 eV. a) and c) (H, K, 0.624)htt plane and 
(H, 0, 0.624)htt scan (circles) for LBCO. b) and d) same for 
LNSCO. In c) and d) the solid line is the result of a fit to a 
Lorentzian squared lineshape. Scattered intensities are nor- 
malized to the respective (002)htt Bragg peaks and therefore 
may be compared directly to each other. 



poorly correlated in LNSCO as compared to LBCO. At 
this point, it is not possible to know which, if any, of 
these explanations is the correct one. 

We next discuss the lattice modulation that is con- 
comitant with the electronic modulation, as measured 
with hard x-rays. In Fig. |3j we show scans performed 
with Ej=11.3 keV through the (0.24,0,11.5) charge stripe 
wavevector. We have again normalized the scattered 
intensity by a nearby Bragg peak, in this case the 
(1, 0, 12)htt- Due to the poorer experimental resolution 
of the hard x-ray measurements, it is necessary to de- 
convolve the experimental resolution. This latter was 
determined from the (1,0, 12) H tt reflection. As an up- 
per bound on the resolution, we took the sharpest of the 
LBCO and LNSCO (1,0,12) H tt reflections. After de- 
convolving the instrumental resolution we find Clb'co = 
178 ± 50 A and ^lnsco = 107 ± 20 A - Comparing the 
widths of the stripe peaks as measured with soft and hard 
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FIG. 3. Comparison of stripe order peaks in a) LBCO and b) 
LNSCO as measured by hard x-rays, E*=11.3 keV at T=12 K. 
Solid lines are the results of fits to Lorentzian squared line- 
shapes. Dashed lines represent the width of the (1, 0, 12)htt 
nearby Bragg reflection in each case, demonstrating the two 
samples are of similar quality. Scattered intensities are nor- 
malized to the (1, 0, 12)htt Bragg peak and therefore may be 
compared directly to each other. 



x-rays, we find that the widths are equal within exper- 
imental error for the two x-ray energies in both LBCO 
and LNSCO. 

From the hard x-ray deconvolved width and the peak 
intensity, we can calculate the ratio of the lattice modula- 
tions for the two systems following the same arguments 
as before. We find that the ratio of the integrated in- 
tensities of the lattice modulations is V T l °f tlce = 20 ± 6, 
i.e. the lattice modulation amplitude is ~ 4.4±0.7 times 
larger in LBCO than in LNSCO. 

Before we turn to a discussion of the significance of 
these results, we present the temperature dependence of 
the respective modulations, and of the LTT phase, for 
the two systems (Fig. Eh. For soft x-rays, the latter is 
measured through the forbidden (OOI)htt reflection, ob- 
servable only in the LTT phase at E t — 532 eV^. For 
hard x-rays, the (1,0, 12)htt reflection, observable only 
in the LTT phase is used. In the case of the superlattice 
modulations, the integrated intensity of the reflection is 
taken to be oc I x T 2 . This is consistent with two di- 



FIG. 4. Temperature dependence of the lattice modulation 
(red circles), electronic modulation (magenta up triangles), 
spin modulation (black crosses) and the LTT phase as mea- 
sured by soft (green down triangles) and hard x-rays (blue 
squares) in both LBCO (a, top panel) and LNSCO (b, bot- 
tom panel). The data presented for the spin modulation were 
taken from Ref . H2l and Ref. OH for LBCO and LNSCO, re- 
spectively. 



mensional stripes, where the L-axis dependence of the 
scattering is dominated by the structure factor due to 
the very short correlation length^!. 

There are several important things to note from this 
plot. First, for LBCO it is clear that an observable peak 
is visible (both the strain modulation and the charge den- 
sity modulation) only within the LTT phase. In addition, 
the integrated intensity is constant, within experimen- 
tal error, for both the electronic and structural modu- 
lations. This demonstrates that the amplitudes of both 
these modulations are constant up to the LTT transition 
temperature, that is the order parameter is not chang- 
ing, only the correlation length. This suggests that had 
the LTT transition occurred at a higher temperature in 
LBCO then it is likely that both the electronic and struc- 
tural modulations would also have persisted to higher 
temperatures. 

In contrast, for LNSCO (Fig. [4|d) the amplitude decays 
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upon warming and is unobservable above Tco=55 K, well 
below the LTT transition at Tltt=68 K. This behavior 
is reminiscent of recent results in Lai.gEuo.2-£cSr :c Cu04 
which has an even larger difference between Tltt and 
T co of ~ 55 kPE3 

Finally, it is worth noting that the temperature de- 
pendence of the LTT-peaks are identical in both samples 
when measured at 11.3 keV and 530 eV, demonstrat- 
ing that the soft x-rays are probing the same bulk-like 
behavior as the hard x-rays and that the samples are 
homogeneous. 

The results presented above represent the first time 
these hole and lattice modulations have been character- 
ized in this detail, and as a result, we are able to draw a 
number of conclusions. First, we address the incommen- 
surability, that is the deviation from the ideal value of 
5 = 0.25 at this doping. In LNSCO this was previously 
attributed to a mixture of two different commensurate 
stripe spacingiPSl. A random ordering of these compo- 
nents gave a good match to the spin and lattice incom- 
mensurabilities and correlation lengths determined with 
neutrons. The considerably longer charge and lattice cor- 
relation lengths in LBCO compared to LNSCO together 
with the same incommensurability in the two systems, is 
not compatible with this random-ordering model. Within 
a model of two integral stripe spacings, correlated order- 
ing of the two components would be necessary to explain 
these results. 

We now address the differences between LBCO and 
LNSCO, searching for hints as to the origin of the differ- 
ences in stripe amplitude. First, the average struc tures 
of LBCO and LNSCO are remarkably similaJ^^'. Sec- 
ondly, the electron-phonon coupling for the Cu-0 bond 
stretching mode appears to be the sam^. Finally, the 
average A-site cation radius is again very similar: 1.232 A 
and 1.211 A for LBCO and LNSCO respectively^. What 
then is the origin of the difference in the two systems? 
Local pinning of stripe domains (as occurs with Zn 
doping can cause neighboring stripe domains to be out 
of phase with one another, and destructive interference 
may limit correlation lengths and average amplitudes. 
The orthorhombic strains (associated with the octahe- 



dral tilts) at temperatures just above the LTT transition 
are twice as large in LNSCO^ as in LBCCff^ and are a 
candidate for the origin of this pinning. Indeed in the or- 
dered phase the LTT tilt angle is larger in LNSCO than 
in LBCO (4.4° compared to 3.5°, see Ref. Hj]) and it 
has previously been shown that the electronic properties 
are sensitive to this degree of tilt, with a critical value 
of $ c = 3.6° below which the electronic properties were 
unaffectecpH By analogy with perovskite manganites^, 
we may expect that the larger modulations of the Cu-O- 
Cu bond angles resulting from this tilt will be associated 
with a greater tendency to localize charge and perturb 
the stripe modulations, especially in the vicinity of diva- 
lent dopant ions. Thus, the reduced average tilt angle in 
LBCO might explain the longer correlation lengths and 
enhanced modulation amplitudes of the periodic stripes. 

In summary, our comparison shows that the in-plane 
correlations of the charge stripes in these two systems 
are isotropic, and a factor of two larger in LBCO. These 
two systems show different temperature dependences, 
with the integrated intensity of the charge modulation 
in LBCO constant up to the LTT transition, while in 
LNSCO it decreases on warming, disappearing at ~ 55 K, 
well below the LTT transition in this material. Compar- 
ing the amplitudes of the lattice modulation we find that 
the lattice modulation in LBCO is 4 times stronger than 
in LNSCO. Performing the same comparison on the elec- 
tronic modulation suggests that the electronic modula- 
tion, as measured by diffraction is 10 times stronger in 
LBCO when compared with LNSCO. This result, based 
on measurement of the integrated intensity of the stripe 
modulation, appears at odds with measurements per- 
formed by resonant lineshape analysis. It is clear that 
further work is required to solve this puzzle, and we hope 
that this work encourages such measurements. 

We gratefully acknowledge D.S. Coburn, W. Leon- 
hardt, W. Schoenig, and S. Wirick for their technical 
support. Work performed at BNL was supported by the 
US Department of Energy, Division of Materials Science, 
under contract No. DE-AC02-98CH10886. J.G. appreci- 
ates the support by the DFG through the Emmy-Noether 
program. 



* swilkins@bnl.gov 

1 J. M. Tranquada, J. D. Axe, N. Ichikawa, et ai, Phys Rev 
Lett 78, 338 (1997). 

2 A. R. Moodenbaugh, Y. Xu, M. Suenaga, et ai, Phys. Rev. 
B 38, 4596 (1988). 

3 M. Fujita, H. Goka, K. Yamada, et al, Phys. Rev. B 70, 
104517 (2004). 

4 S. Kivelson, E. Fradkin, and V. Emery, Nature 393, 550 
(1998). 

5 V. J. Emery, S. A. Kivelson, and O. Zachar, Phys. Rev. 
B 56, 6120 (1997). 

6 Q. Li, M. thicker, G. D. Gu, et al, Phys Rev Lett 99, 
067001 (2007). 

7 J. M. Tranquada, G. D. Gu, M. Hiicker, et al, Phys. Rev. 



B 78, 174529 (2008). 

8 M. J. Lawler, K. Fujita, J. Lee, et al, Nature 466, 347 
(2010). 

9 R. Daou, J. Chang, D. Leboeuf, et ai, Nature 463, 519 
(2010). 

10 P. Abbamonte, A. Rusydi, S. Smadici, et ai, Nature Phys. 
1, 155 (2005). 

11 Y.-J. Kim, G. D. Gu, T. Gog, et al, Phys. Rev. B 77, 
064520 (2008). 

12 M. Hiicker, M. von Zimmermann, G. Gu, et al, Phys. Rev. 
B 83, 104506 (2011). 

13 J. Tranquada, B. Sternlieb, J. D. Axe, et ai, Nature 375, 
561 (1995). 

14 J. Fink, E. Schierle, E. Weschke, et ai, Phys. Rev. B 79, 



6 



100502 (2009). 

15 J. Fink, V. Soltwisch, J. Geek, et al, Phys. Rev. B 83, 
092503 (2011). 

16 C. T. Chen, F. Sette, Y. Ma, et al, Phys Rev Lett 66, 104 
(1991). 

17 H. Romberg, M. Alexander, N. Nucker, et al, Phys. Rev. 
B 42, 8768 (1990). 

18 S. Partzsch and S. B. Wilkins, Unpublished. 

19 J. M. Tranquada, J. D. Axe, N. Ichikawa, et al, Phys. Rev. 
B 54, 7489 (1996). 

20 A. Rusydi, P. Abbamonte, H. Eisaki, et al., Phys Rev Lett 
97, 016403 (2006). 

21 L. M. Rodriguez-Martinez and J. P. Attfield, Phys. Rev. 
B 54, R15622 (1996). 

22 B. Nachumi, Y. Fudamoto, A. Keren, et al., Phys. Rev. B 
58, 8760 (1998). 

23 M. von Zimmermann, A. Vigliantc, T. Niemoller, et al., 



Europhys. Lett. 41, 629 (1998). 

24 J. M. Tranquada, N. Ichikawa, and S. Uchida, Phys. Rev. 
B 59, 14712 (1999). 

25 J. D. Axe and M. K. Crawford, Journal of Low Tempera- 
ture Physics 95, 271 (1994). 

26 S. Katano, J. Fernandez-Baca, S. Funahashi, et al., Physica 
C: Superconductivity 214, 64 (1993). 

27 D. Reznik, L. Pintschovius, M. Ito, et al, Nature 440, 1170 
(2006). 

28 R. D. Shannon, Acta Crystallographica Section A 32, 751 
(1976). 

29 H. Kimura, M. Kofu, Y. Matsumoto, et al, Phys Rev Lett 
91, 067002 (2003). 

30 B. Buchner, M. Breuer, A. Freimuth, et al, Phys Rev Lett 
73, 1841 (1994). 

31 P. G. Radaelli, M. Marezio, H. Y. Hwang, et al, Phys. 
Rev. B 54, 8992 (1996). 



